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Abstract 

OS ' 

vVe present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneus reionization of the universe 
based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at I = 2000 — 8000 
£ — with maximum values of [£(£ + l)CV/(2-7r)] max ~ 4 — 7 x 10 -13 . The scale roughly corresponds to the typical ionized bubble sizes 
observed in our simulations, of ~ 5 — 20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal 
above t = 3000. At large scales the patchy kSZ signal depends only on the source efficiencies, ft is higher when sources are more 
efficient at producing ionizing photons, since such sources produce larger ionized regions, on average, than less efficient sources. The 

(— I introduction of sub-grid gas clumping in the radiative transfer simulations produce significantly more power at small scales, but 

O fias little effect at large scales. The patchy reionization kSZ signal is dominated by the post-reionization signal from fully-ionized 
1 gas, but the two contributions are of similar order at scales I ~ 3000 — 10 4 , indicating that the kSZ anisotropies from reionization 

Q are an important component of the total kSZ signal at these scales. 
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5t Introduction 

The secondary anisotropies of the Cosmic Microwave 
Background (CMB) are emerging as one of the most pow- 
erful tools in cosmology. The small-scale anisotropies are 
probes of cosmological structures and are thus a valuable 
tool for studying their formation and properties. The key 
effect resulting in small-scale CMB anisotropies is the 
Sunyaev-Zel'dovich (SZ) effect (Zeldovich and Sunyaev, 
1969), produced by Compton scattering of the CMB pho- 
tons on moving free electrons. When this is due to ther- 
mal motions the effect is referred to as thermal SZ effect 
(tSZ), while when it is due to the electrons moving with 
a net bulk peculiar velocity it is called kinetic SZ effect 
(kSZ) (Sunyaev and Zeldovich, 1980). The former, com- 
ing largely from hot gas in low-redshift galaxy clusters, 
is the dominant effect, but has a different spectrum than 
the CMB primary anisotropies, with a characteristic zero 
at ~ 217 MHz, so its contribution can be separated. By 
contrast, kSZ anisotropies have a spectrum identical with 
that of the primary anisotropies. Separation of its effects 



are possible only by virtue of its differing spatial struc- 
ture. We consider the kSZ fluctuations as composed of 
two basic contributions, one coming from inhomogeneous 
reionization and the other from the fully-ionized gas after 
reionization. The calculation of the patchy reionization 
contribution to the kSZ signal based on the first large-scale 
radiative transfer simulations of reionization is the main 
focus of this paper. This signal provides a signature of the 
character of reionization that will complement other ap- 
proaches like observations of the redshifted 21-cm line of 
hydrogen and surveys of high-z Ly-a emitters (e.g. Tozzi 
et al, 2000; Iliev et al, 2002; Rhoads et al., 2003; Stanway 
and et al., 2004; Gnedin and Prada, 2004; Furlanetto et al., 
2004a; Mellema et al., 2006b; Malhotra and Rhoads, 2005; 
Shapiro et al., 2006; Bunker et al., 2006). 

The kSZ from a fully-ionized medium has been studied 
previously by both analytical and numerical means. When 
linear perturbation theory is used, the effect is associated 
with the quadratic nonlinearities in the electron density 
current and is usually referred to as the Ostriker-Vishniac 
effect (Ostriker and Vishniac, 1986; Vishniac, 1987). Cal- 



culating the full, nonlinear effect is much more difficult to 
do analytically, although some models have been proposed 
(e.g. Jaffe and Kamionkowski, 1998; Hu, 2000; Ma and 
Fry, 2002; Zhang et al., 2004). This effect involves coupling 
of very large-scale to quite small-scale density fluctuation 
modes and thus requires a very large dynamic range in order 
to be simulated correctly. Current simulations are in rough 
agreement, but have not quite converged yet (e.g. Springcl 
et al., 2001; Ma and Fry, 2002; Zhang et al., 2004). All sim- 
ulations predict significant enhancement of the small-scale 
anisotropies compared to the analytical theory, due to the 
nonlinear evolution. 

Calculations of the second important kSZ contribution, 
from patchy reionization, are even more complicated. In 
addition to the above difficulties, it also requires detailed 
modelling of the radiative transfer to derive the sizes and 
distributions of H II regions in space and time (the reion- 
ization geometry) and how these correlate with the veloc- 
ity and density fields. To date this problem has been stud- 
ied by only a few recent works. Most of these estimates 
were done by semi-analytical models, which used different 
simplified approaches to model the inhomogeneous reion- 
ization (Gruzinov and Hu, 1998; Santos et al., 2003; Zahn 
et al., 2005; McQuinn et al., 2005). While such models are 
useful since they are simple and much cheaper to calculate 
than full radiative transfer simulations, and are thus con- 
venient for exploring various scenarios, their results must 
be checked against full detailed simulations in order to as- 
certain their reliability. The only two existing numerical 
studies of this effect, Gnedin and Jaffe (2001) and Sal- 
vaterra et al. (2005) used fairly small computational boxes 
(4h _1 Mpc and 20h~ 1 Mpc, respectively), and, as a conse- 
quence, the results they derived significantly underestimate 
the signal since such small boxes do not include the larger- 
scale velocity modes which are quite important. 

There are several upcoming experiments which would 
search for kSZ effect, in particular the Atacama Cosmol- 
ogy Telescope (ACT) 1 and South Pole Telescope (SPT) 2 , 
both of which are expected to be operational by early 2007. 
These experiments would have ~ 1' resolution and ~ /xK 
sensitivities, which should be sufficient to detect the kSZ 
effect. 

In this paper we present the first calculations of the kSZ 
effect from patchy reionization based on large-scale radia- 
tive transfer simulations. We simulate reionization using 
a (lOOh^Mpc) 3 simulation volume, which is sufficient to 
capture the relevant large-scale density and velocity pertur- 
bations, an important improvement over previous efforts. 
We use the ray-tracing code C 2 — Ray (Mellema et al., 
2006a) to follow the radiation from all ionizing sources in 
that volume identified with the resolved halos, which are of 
dwarf galaxy size or larger. The halos and underlying den- 
sity field are provided by a very large N-body simulation 
with the code PMFAST (Merz et al., 2005). Our methodol- 
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ogy, simulation parameters and results on the reionization 
character, geometry and observability of the rcdshiftcd 21- 
cm line of hydrogen were discussed in detail in Iliev et al. 
(2006, Paper I) and Mellema et al. (2006b, Paper II). 

2. Simulations 

Our simulations follow the evolution of a comoving sim- 
ulation volume of (100 /i _1 Mpc) 3 , corresponding to an an- 
gular size <~ 1 deg. Our radiative transfer grids are of sizes 
203 3 or 406 3 , which allows us to derive reliably the angular 
sky power spectra for/ w 430-90, 000 (180, 000). We utilize 
the flat-sky approximation, which is appropriate for such 
relatively small angular sizes. Our simulation methodology 
and parameters were described in detail in Papers I and II. 
Here we provide just a brief summary. We start with per- 
forming a very large pure dark matter simulation of early 
structure formation, with 1624 3 particles and 3248 3 grid 
cells using the code PMFAST (Merz et al., 2005). This al- 
lows us to reliably identify (with 100 particles or more per 
halo) all halos with masses 2.5 x 10 9 M Q or larger. We find 
and save the halo catalogues, which contain the halo po- 
sitions, masses and detailed properties, in up to 100 time 
slices starting from high redshift (z ~ 30) to the observed 
end of reionization at z ~ 6. We also save the correspond- 
ing density and bulk peculiar velocity fields at the resolu- 
tion of the radiative transfer grid. Unfortunately, radiative 
transfer simulations at the full grid size of our N-body com- 
putations are impractical on current computer hardware. 

All identified halos are assumed to be sources of ionizing 
radiation and each is assigned a photon emissivity propor- 
tional to its mass, M, according to 

N~ = f~ M \ , (1) 
fim p t s tt 

where t s is the source lifetime, m p is the proton mass, and / 7 
is a photon production efficiency which includes the num- 
ber of photons produced per stellar atom, the star forma- 
tion efficiency (i.e. what fraction of the baryons are con- 
verted into stars) and the escape fraction (i.e. how many 
of the produced ionizing photons escape the halos and are 
available to ionize the IGM) . 

The radiative transfer is followed using our fast and 
accurate ray-tracing photoionization and non-equilibrium 
chemistry code C 2 — Ray. The code has been tested in 
detail for correctness and accuracy against available an- 
alytical solutions and a number of other cosmological 
radiative transfer codes (Mellema et al., 2006a; Iliev and 
et al., 2006). The radiation is traced from every source on 
the grid to every cell. 

We have performed four radiative transfer simulations. 
These share the source lists and density fields given by the 
underlying N-body simulation, but adopt different assump- 
tions about the source efficiencies and the sub-grid density 
fluctuations. The runs and notation are the same as in Pa- 
per II: runs f2000 and f250 assume / 7 = 2000 and 250, 



respectively, and no sub-grid gas clumping, while f2000C 
and f250C adopt the same respective efficiencies, / 7 = 2000 
and 250, but also add a sub-grid gas clumping, C(z) = 
(n 2 )/{n) 2 , which evolves with redshift according to 

C subgri dW - 27.466 e -°- 114z + - 00132822 . (2) 

The last fit was obtained from another high-resolution PM- 
FAST N-body simulation, with box size (3.5 h -1 Mpc) 3 and 
a computational mesh and number of particles of 3248 3 and 
1624 3 , respectively. These parameters correspond to a par- 
ticle mass of 10 3 Mq and minimum resolved halo mass of 
10 5 M Q . This box size was chosen so as to resolve the scales 
most relevant to the gas clumping - on scales smaller than 
these the gas fluctuations would be below the Jeans scale, 
while on larger scales the density fluctuations are already 
present in our computational density fields and should not 
be doubly-counted. The expression in equation (2) excludes 
the matter inside collapsed minihalos (halos which are too 
small to cool atomicly, and thus have inefficient star forma- 
tion) since these are shielded, unlike the generally optically- 
thin IGM. This self-shielding results in a lower contribu- 
tion of the minihalos to the total number of recombinations 
than one would infer from a simple gas clumping argument 
(Shapiro et al., 2004; Iliev ct al., 2005). The effect of mini- 
halos could be included as sub-grid physics as well, see Cia- 
rdi et al. (2006). This results in slower propagation of the 
ionization fronts and further delay of the final overlap. The 
halos that can cool atomicly are assumed here to be ioniz- 
ing sources and their recombinations are thus implicitly in- 
cluded in the photon production efficiency / 7 through the 
corresponding escape fraction. 

3. kSZ from Patchy Reionization 

The kSZ effect is the CMB temperature anisotropy along 
a line-of-sight (LOS) defined by a unit vector n induced by 
Thomson scattering from flowing electrons: 

AT f 

- / drje- T ^an e a T n ■ v, (3) 

-L CMB J 

where 77 = J Q dt'/a(t') is the conformal time, a is the scale 
factor, <jt — 6.65 x 10~ 25 cm~ 2 is the Thomson scattering 
cross-section, and r is the corresponding optical depth. 

We calculate the kSZ anisotropy signal from our simu- 
lation data as follows. We first calculate the line-of-sight 
integral in equation (3) for each individual output time- 
slice of the radiative transfer simulation and for all LOS 
along each of the three box axes. The contribution to the 
total LOS integral from each light-crossing time of the box 
is then obtained by linear interpolation between the near- 
est results from the simulation output times. The individ- 
ual light-crossing time contributions are then all added to- 
gether to obtain the full LOS integral given in equation (3). 
All these integrals are done for each LOS through the box 
along the direction of light propagation, which allows us 
to produce kSZ maps, in addition to the statistical signals. 



In order to avoid artificial amplification of the fluctuations 
resulting from repeating the same structures along the line 
of sight, after each light-crossing time we randomly shift 
the box in the directions perpendicular to the LOS and by 
rotating the box, so that the LOS cycles x, y and z axes of 
the simulation volume. 

In addition to our simulations, for comparison we also 
consider two simplified cases, of instant and of uniform 
reionization. We define instant reionization as a sharp tran- 
sition from completely neutral to fully-ionized IGM at red- 
shift instant- We pick Zi ns t a nt = 13, which yields the same 
integrated electron scattering optical depth as our simula- 
tion f250. We also define a "uniform reionization" scenario 
to be one that has the same time-dependent reionization 
history (and hence, also the same r es ) as simulation f250, 
but spatially uniform (i.e. not patchy). We then use the 
same density and velocity data and same procedures as for 
the actual simulations to derive the kSZ temperature fluc- 
tuations for these two scenarios. We consider these simple 
models in order to demonstrate the effects of reionization 
being extended and patchy in nature. 

4. Results 

We show the kSZ maps of temperature fluctuations, 
ST/Tqmb, yielded by all of our cases in Figures 1 (f2000 
and f250), 2 (f2000C and f250C), and 3 (instant and uni- 
form reionization). These have a total angular size of ap- 
proximately 50' x 50', corresponding to our computational 
box size. The resolution of the maps is that of the full 
radiative transfer grid (203 x 203 pixels), corresponding 
to pixel resolution of ~ 0.25'. All maps utilize the same 
color map in order to facilitate their direct comparison. 
The maps deriving from our simulations of inhomogeneous 
reionization (Figures 1 and 2) all show fairly strong fluctu- 
ations, both positive and negative, of order ST ~ 10 \iK at 
angular scales of a few arcminutes (~ 10 h _1 Mpc) and up 
to ST ~ 20 \iK at the full map resolution. The fluctuations 
are at somewhat smaller scales when the sources are less 
efficient photon producers (f250 and f250C), compared to 
the high-efficiency cases (f2000 and f2000C). The varia- 
tions are noticeably enhanced when sub-grid gas clumping 
is included in the reionization model (f2000C and f250C) 
and there are a number of regions with very strong fea- 
tures. In comparison, the artificial models of instant and 
uniform reionization (Figure 3) show fluctuations with 
much lower amplitudes and with a less well-defined typi- 
cal scale. Since these simple scenarios were constructed to 
produce the same total electron scattering optical depth 
as simulation f250 and they use the same density and ve- 
locity fields as the simulations, any differences we see are 
due to the patchiness of reionization. The lower signal in 
these last two cases is expected since the kSZ effect in 
uniformly-ionized gas exactly cancels in the first order. 
This cancellation is broken when patchiness is present, 
enhancing the anisotropies. 




Fig. 1. kSZ maps from simulations f2000 (left) and f250 (right). (Images produced using the Ifrit visualization package of N. Gnedin). 
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Fig. 2. kSZ maps from simulations f2000C (left) and f250C (right). 
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Fig. 3. kSZ maps from simulations: (left) assuming instant reionization at Zinstant = 13 which gives the same total electron scattering optical 
depth as simulation f250, and (right) spatially-uniform reionization with the same reionization history and thus same total electron scattering 
optical depth as simulation f250. 



The mean and rms of the temperature fluctuations for 
all cases are summarized in Table 1. Based on these, we 
see that, indeed, the uniform and instant reionization cases 
have means very close to zero, while the realistic, patchy 
reionization simulations yield mean temperature fluctua- 
tions that are low, of order ~ 10~ 7 , but not zero. The rms 
values are ~ 10~ 6 for all patchy reionization cases and lower 



than that by factor of ~ 2 (~ 3) for the instant (uniform) 
reionization cases. 

These observations based on the maps are confirmed by 
their corresponding pixel PDF distributions, shown in Fig- 
ures 4 and 5. On each panel we also plot the Gaussian distri- 
bution with the same mean and standard deviation. All dis- 
tributions are surprisingly close to Gaussian around their 
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Table 1 

Mean and rms values for 5Tksz /^CMB • 

simulation f2000 f2000C f250 f250C uniform instant 

mean 8.86 x 1CT 8 9.88 x 1CT 8 -2.17 x 10~ 7 5.57 x 10" 8 -3.29 x 10" 12 -3.37 x 10" 12 

rms 1.13 x 10~ 6 1.22 x 10~ 6 1.06 x 10~ 6 1.09 x 1CT 6 4.24 x 10~ 7 6.67 x 10~ 7 
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Fig. 4. PDF distribution of <57ksz/7cMB (solid) vs. Gaussian distri- 
bution with the same mean and width (dotted). 

mean values given the maps. However noticeable depar- 
tures from Gaussianity do occur in the wings of the PDFs. 
In particular, the PDFs for patchy reionization with sub- 
grid clumping (f2000C, and f250C) are significantly non- 
Gaussian. For the realization simulated there is an over- 
abundance of bright regions by up to an order of mag- 
nitude compared to the corresponding Gaussian distribu- 
tions. The reionization scenarios without sub-grid clump- 
ing show much weaker non-Gaussian features there. This 
indicates that the observed PDF of kSZ from patchy reion- 
ization may give us important information on the level of 
small-scale gas dumpiness during reionization. The PDFs 
derived from the uniform and instant reionization scenarios 
are significantly less wide than the simulated ones (which 
was also shown by their lower rms values, as was discussed 
above) , and are also much closer to Gaussian. 

The kSZ anisotropy signal from inhomogeneous reion- 
ization dominates the primary CMB anisotropy above £ = 
3000. The sky power spectra (Figure 6) peak strongly at 
£ = 3000 - 5000, to a maximum of [£(£ + l)Q/(27r)] max - 
7 x 10 -13 when the ionizing sources are highly-efficient 
(cases f2000 and f2000C). The peak values for the simu- 
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Fig. 5. PDF distribution of ST^sz/Tcmb (solid) vs. Gaussian distri- 
bution with the same mean and width (dotted). 

lations with lower-efficiency sources (f250 and f250C) are 
only slightly lower, at [£{£ + l)Ct/(2w)] max - 4x 10~ 13 , 
and the peaks are somewhat broader and moved to smaller 
scales, £ ~ 3000 — 7000, than for the high-efficiency cases. 
The introduction of sub-grid gas clumping in the radiative 
transfer simulations produces significantly more power at 
small scales and slightly broader peaks when compared to 
the similar cases without sub-grid clumping. The scales at 
which the power spectra peak roughly correspond to the 
typical ionized bubble sizes observed in our simulations, 
namely ^5 — 20 Mpc comoving. These typical sizes depend 
on the assumed source efficiencies and sub-grid clumping. 
At large scales the patchy kSZ signal depends only on the 
source efficiencies, being higher for more efficient photon 
emitters, since these tend to produce larger ionized regions 
on average. 

In contrast, the uniform reionization scenario (with the 
mean reionization history of simulation f250) yields a kSZ 
signature which is much lower than the non-uniform reion- 
ization scenarios, indicating a very large boost of the signal 
due to the patchiness. The boost is largest, approximately 
one order of magnitude, at and above the typical scale of 
the patches (£ < 7000), but is still exists at smaller scales, 
where it is a factor of two or more. The other simplified 
scenario, of instant reionization with the same integrated 
optical depth at f250, produces larger kSZ anisotropy than 
a uniform reionization does. However, it is still well below 
the realistic patchy reionization signals, by factor of ~ 3 for 
£ < 7000. At smaller scales the reionization signal from the 
instant reionization scenario becomes similar to the ones 
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Fig. 7. Contributions to the total kSZ signal (blue) from different redshifts for case f250 (left), uniform (middle) and instant reionization 
(right). Shown are the signals for every three box light crossing times, roughly corresponding to (bottom to top on the left) redshifts z > 15 
(green), 15 > z > 11 (red), z < 11 (magenta). In the last two cases the z > 15 contribution is very low (uniform) or zero (instant) and thus 
not show. 
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Fig. 6. Sky power spectra of <57kgz/TcMB fluctuations resulting 
from our simulations: f2000 (black), f250 (blue), f2000C (red) and 
f250C (green). For comparison, we also show the results from simple 
models which utilize the same density and velocity fields as the actual 
simulations, but different assumptions about the gas ionization in 
space and time, a uniform reionization with the same reionization 
history as simulation f250 (dark red, short-dashed) and an instant 
reionization model with the same integrated optical depth T ea as 
simulation f250 (cyan, short-dashed). The primary CMB anisotropy 
signal is also shown (brown, dotted). 

from simulated non-uniform reionization. The distribution 
of the power in these two simplified scenarios is much flat- 
ter, with less indication of a characteristic scale. This im- 
plies that the sharp peaks yielded by the inhomogeneous 
reionization scenarios are dictated by the size of the ion- 
ized patches, since both the density and the velocity fields 



are shared among the models. 

This behaviour can be understood further by consider- 
ing the contributions from different redshift intervals to the 
integrated signal. In Figure 7 we show the contributions to 
the total signal from different redshift intervals for simula- 
tion f250, as well as the uniform and instant reionization 
scenarios. We plot the contribution from every three light- 
crossing times of the box, corresponding roughly to z > 
15, 15 > £ > 11 and z < 11. The mass- weighted ionized 
fraction in the highest redshift interval is x m < 0.01, and 
consequently its contribution to the kSZ effect (case f250, 
bottom curve) is low, with a maximum of ~ 10~ 14 . This 
high-redshift contribution is strongly peaked at very small 
scales (£ > 10 4 ), reflecting the fact that at that time the 
typical ionized bubble is still quite small, less than few Mpc 
in size. The contribution from the middle redshift interval, 
11 < z < 15, (corresponding to 0.8 > x m > 0.01) peaks 
at roughly the same scales as the integrated signal. It con- 
tributes about half of the total signal around the peak, but 
a smaller fraction at scales above and below that. A half 
or more of the integrated kSZ signal at any scale is con- 
tributed by the lowest redshifts (z < 11, x m ^ 0.8). The 
low- redshift power spectrum is fairly flat, with a weak peak 
at relatively large scales, I ~ 2000. 

In the homogeneous (uniform and instant) reionization 
cases the shape of the power spectra is essentially the same 
for all redshift intervals, quite flat, with a broad peak at 
I ~ 5000. The signal from the uniform reionization scenario 
is completely dominated by the lowest redshifts, with less 
than a few per cent coming from the middle redshift in- 
terval, and essentially no contribution from high redshifts. 
Since the velocity and density fields are shared, the instant 
reionization scenario contributions differ from the uniform 
scenario ones only by being weighted by the ionized frac- 
tion in the uniform case. The instant reionization scenario 
assumes full ionization after instant = 13 and fully-neutral 
gas before that. Thus, as a consequence to the kSZ signal 
being stronger at later times, the instant reionization inte- 



grated signal is higher, as well. 
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Fig. 8. Sky power spectra of ST^sz/Tqmb fluctuations from our 
simulations f2000 (black), f250 (blue), f2000C (red) and f250C 
(green) compared to the after-reionization kSZ signals (assuming 
overlap at z ov = 8): linear Ostriker-Vishniac effect, labelled 'OV 
(long-dashed, pink), the same, but using the nonlinear power spec- 
trum of density fluctuations, labelled 'NL' (long-dashed, dark red) 
and a fully-nonlinear model matched to high-resolution hydrody- 
namic simulations of Zhang et al. (2004), labelled 'Zhang et al' 
(short-dashed, dark green). The primary CMB anisotropy signal is 
also shown (brown, dotted). 

Finally, we compare the two contributions to the total 
kSZ anisotropies, the patchy reionization signal we have 
calculated and the kSZ signal from fully-ionized gas after 
reionization. In Figure 8 we show our predictions against 
three representative calculations of the post-reionization 
contribution: the quadratic-order Ostriker-Vishniac effect 
(Vishniac, 1987) expressed in terms of a product of the lin- 
ear power spectrum; the same expression, but with the non- 
linear density power spectra substituted for one of the lin- 
ear ones, which partially accounts for the nonlinear effects; 
and the recent detailed nonlinear model of Zhang et al. 
(2004). The last three models are rescaled to our adopted 
cosmology using £(£ + 1)CV oc erf (Zhang et al., 2004), and 
assume reionization overlap at z ov = 8, to allow for direct 
comparison. Confronting first the three uniform-ionization 
models, we see that they agree fairly well at large, linear 
scales, but strongly diverge, by up to one order of magni- 
tude at small scales, where the nonlinearities become very 
important. Including the corrections due to the nonlinear 
density power spectrum yields modestly higher kSZ signal 
than the linear OV calculation, and similar power spec- 
trum shape overall. Both peak at roughly the same scale, 
£ ~ 2000. In contrast, the (Zhang et al., 2004) result finds 



much larger signal, especially at small scales, which peaks 
at £ ~ 20, 000. Compared to our patchy reionization kSZ 
predictions, all uniform-ionization power spectra are much 
less sharply peaked, i.e. they lack a well-established charac- 
teristic scale. In terms of their peak values, the OV and OV 
with nonlinear corrections models yield values similar to 
our patchy reionization results, while (Zhang et al., 2004) 
find almost a factor of 2 higher signal at the patchy sig- 
nal peak, and even more pronounced dominance at smaller 
scales. 

5. Summary and Discussion 

We have derived the kSZ CMB anisotropies due to the 
inhomogeneous reionization of the universe. This is the 
first such calculation based on detailed, large-scale radia- 
tive transfer simulations of this epoch. The resulting sky 
power spectra peak at I = 2000 — 8000 with maximum val- 
ues of [£(l + 1 )Ce/(2n)] maK - 4 - 7 x 10~ 13 . This roughly 
corresponds to the typical ionized bubble sizes observed 
in our simulations, which is ^ 5 — 20 Mpc, depending on 
the assumed source efficiencies and the gas clumping at 
very small scales. The kSZ anisotropy signal from reioniza- 
tion dominates the primary CMB signal above £ — 3000. 
At large scales the patchy kSZ signal depends only on the 
source efficiencies. It is higher when sources are more effi- 
cient at producing ionizing photons, since such sources pro- 
duce large ionized regions, on average, than less efficient 
sources. The introduction of sub-grid gas clumping in the 
radiative transfer simulations produce significantly more 
power at small scales, but has little effect at large scales. 
The integrated kSZ signal is strong enough to be detected 
by upcoming experiments, like ACT and SPT, but it seems 
difficult to separate the patchy reionization signal from the 
contribution by the fully-ionized gas after reionization. 

Previous simulations of the kSZ effect from patchy reion- 
ization (Gnedin and Jaffe, 2001; Salvaterra et al., 2005) 
predicted significantly lower kSZ signals than the ones we 
found. Their power spectra reach maximum values of [£(£ + 
l)Q/(27r)] max - 2 x 10~ 14 and - 1.6 x 10~ 13 , respectively, 
compared to [£(£ + l)Ci/(2n)] max - 4 - 7 x 10~ 13 for our 
simulations. This discrepancy is due to the small volumes 
used in these simulations, which significantly reduces the 
power in both the density and the velocity field fluctua- 
tions. At the scale of the simulation box the bulk veloci- 
ties are zero by definition, and the density is the mean one 
for the universe, thus any larger-scale fluctuations are not 
included, and the ones close to the box size are underesti- 
mated. In terms of the calculation of kSZ anisotropy from 
patchy reionization, this problem is more severe for the ve- 
locity field, since its power grows towards large scales, un- 
like the density field fluctuations, which instead decrease at 
large scales. The reionization patchiness imprints its char- 
acteristic scale on the density and velocity fluctuations, ef- 
fectively smoothing the small-scale fluctuations below the 
typical bubble sizes. On the other hand, the large-scale flue- 



tuations still should contribute to the kSZ signal, since the 
H II regions are moving with the large scale bulk motions. 

In addition to the missing large-scale power, these earlier 
reionization simulations also did not follow sufficient vol- 
umes to properly sample the size distribution of the ionized 
bubbles. The typical sizes of these ionized patches are of or- 
der 5-20 comoving Mpc, and even larger at later times, and 
hence require simulation volumes of order (100 h _1 Mpc) 3 
for correct sampling. As a result, the kSZ power spectra 
found by the smaller-box simulations were largely flat, with 
no clear characteristic scale, since they do not reach the 
large characteristic bubble scales, which are dictated by the 
strong clustering of ionizing sources at high redshifts. 

Several semi-analytical models for calculating the kSZ 
signature of patchy reionization have been proposed in 
recent years (Gruzinov and Hu, 1998; Santos et al., 2003; 
McQuinn et al., 2005). Gruzinov and Hu (1998) proposed 
a very simple model, whereby the ionized patches are ran- 
domly distributed and have a given characteristic size R. 
Furthermore, they assumed that the density, velocity and 
ionization fraction fields are all uncorrelated with each 
other. The ionized fraction auto-correlation function is 
approximated as a Gaussian with rms given by the (fixed) 
characteristic H II region size. Under these assumptions 
the kSZ power spectrum can be calculated analytically. 
The quantitative details of the signal predicted by this 
model depend on the assumptions made about the typical 
ionized patch size, and the time and duration of reioniza- 
tion, but generically it predicts a distribution which is very 
strongly peaked, more so than our simulations find. This 
is related to the assumed Gaussian distribution around a 
fixed characteristic size of the ionized patches. The actual 
simulations also find a characteristic size for the ionized 
bubbles, but one that is evolving with redshift, and differ- 
ent size distributions around it, resulting in somewhat less 
sharp peaks than this analytical model predicts. 

Santos et al. (2003) proposed another simple model for 
the reionization patchiness, which assumes that the corre- 
lation function of the ionized patches is proportional to the 
one for the density field, with a bias factor which is time- 
dependent, but not scale-dependent. The power is filtered 
at small scales, below the typical size of the ionized bub- 
bles using a Gaussian filer in k-space. The time-dependence 
of the typical size of the ionized patches is assumed to be 
proportional to (1 — x m (t)) _1/ ' 3 , where x m (t) is the mean 
mass- weighted ionized fraction at time t. The last function 
is based on the semi-analytical models of reionization of 
Haiman and Holder (2003). The resulting kSZ power spec- 
tra are fairly flat, with much less well-defined characteris- 
tic scale than our simulation results, as should be expected 
based on the quickly-evolving typical patch size assumed 
in this model. The peak amplitudes of the results is some- 
what higher than the ones we find, by factors of ~ 2 — 3, 
probably due to their assumed values of the effective bias. 

More recently, McQuinn et al. (2005) applied the semi- 
analytical reionization model of Furlanetto et al. (2004b) 
to estimate the kSZ signal. They found lower fluctuations 



in their "single reionization episode" models (by factor 
of ~ 2), although they have similar mean reionization 
histories (i.e. evolution of the mean ionized fraction) to 
our simulations. Their results predict a peak at somewhat 
larger scales (£ -2000) than ours (£ -2000-9000). Their 
"extended reionization" scenarios peak higher (in some 
cases slightly above [£(£ + l)CV/(27r)] max — 10~ 12 ) again 
at £ —2000. We note that these extended scenarios have a 
very different assumed reionization histories than any of 
our simulations. They found that the kSZ post-reionization 
contribution dominates the patchy signal at the scales of 
interest, in rough agreement with our results. 
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